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A method to calculate two-loop self-energy diagrams of the Standard Model is demon- 
strated. A direct physical application is the calculation of the two-loop electroweak 
contribution to the anomalous magnetic moment of the muon \{g — 2) . Presently we 
confine ourselves to a "toy" model with only fi, 7 and a scalar particle (Higgs). The 
algorithm is implemented as a package of computer programs in FORM. For generating 
and automatically evaluating any number of two-loop self-energy diagrams, a special C- 
program has been written. This program creates the initial FORM-expression for every 
diagram generated by QGRAF, executes the corresponding subroutines and sums up the 
final results. 

Keywords: Standard Model, Feynman diagram, recurrence relations, anomalous mag- 
netic moment 

1. Introduction 

Recent high precision experiments to verify the Standard Model of electroweak 
interactions require, on the side of the theory, high precision calculations resulting 
in the evaluation of higher loop diagrams. For specific processes thousands of multi- 
loop Feynman diagrams do contribute, and it turns out impossible to perform these 
calculations by hand. That makes the request for automatization a high priority 
task. In this direction, several program packages are elaborated (0 and reports 
of the Minami-Tateya collaboration, Bauberger et al., Pukhov et al., Vermaseren 
and J.X.Wang on this conference). It appears absolutely necessary that various 
groups produce their own solutions of handling this problem: the various ways will 
be of different efficiency, have different domains of applicability, and last but not 
least, should eventually allow for completely independent checks of the final results. 
This point of view motivated us to seek our own way for automatic evaluation of 
Feynman diagrams. We have in mind only higher loop calculations (no multipoint 
functions). 

We demonstrate here the functioning of a C-program (TLAMM) for the evalu- 
ation of the two-loop anomalous magnetic moment (AMM) of the muon h{g — 2) , 
but on the same lines we have developed the program to an extent that it can be ap- 
plied to arbitrary processes in the Standard Model. This piloting C-program must 
read the diagrams generated by QGRAF LI for a given physical process, generate 
the FORM □ source code, start the FORM interpreter, read and sum up the results 
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for the class of diagrams under consideration. For the purpose of demonstration, 
here we apply TLAMM to a closed subclass of diagrams of the Standard Model, 
which we refer to as a "toy" model. 

Recent papers reduced the theoretical uncertainty of the muon AMM by partially 
calculating the two-loop electroweak contributions El, B. Their results have been 
obtained in the following approximation: terms suppressed by (1 — 4 sin 2 0) were 
omitted; the fermion masses of the first two families are set to zero; in the third 
family the r- and 6-quark masses are taken to be zero as well; diagrams with two or 

m 2 

more scalar couplings to the muon are suppressed by the ratio jfe and have been 
discarded; the Kobayashi-Maskawa matrix is assumed to be unity; the mass of the 
Higgs particle is large compared to Mz,w- 

All of these approximations, except possibly the last one, are well justified and 
give rise to small corrections only. We consider it of great interest to study also the 
case Mh ~ Mz.w- To perform this calculation is our main physical motivation. 
Apart from that, for technical reasons, it may be interesting to study the functioning 
of TLAMM by calculating all 1832 two-loop diagrams without any approximation. 

The calculation of the anomalous magnetic moment of the muon reduces, after 
differentiation and contractions with projection operators, to diagrams of propaga- 
tor type with external momentum on the muon mass-shell (for details see 0). 

The method of Taylor expansion of diagrams in external momenta and Pade 
approximation u yields only slowly converging results at the threshold. Since in the 
case of the muon AMM we are confronted with a low energy problem it appears 
natural to expand with respect to the heavy masses (Mz, Mh and mt) of the theory. 
The applicability of the "asymptotic-R-operation" i in the limit of large masses 
has to be investigated for diagrams evaluated on the muon mass-shell(i.e. p 2 = 
—mV). Some diagrams already had a threshold at the muon mass-shell before 
the expansion. In other diagrams this threshold appears in some terms of the 
expansion. In dimensional regularization, threshold singularities (like any other 
infrared singularities if they are strong enough) manifest themselves as poles in e 
(d=4-2e). They ought to cancel for the total AMM. We check this in our toy model. 

2. Large-mass expansion 

The asymptotic expansion in the limit of large masses is defined as 

F G (q,M,m,e) M ~°° £ F Gh (q, m, e) o T^^F^q* , M, m\ e) (1) 

7 

where G is the original graph, 7's are subgraphs involved in the asymptotic expan- 
sion, G/7 denotes shrinking 7 to a point; Fy is the Feynman integral corresponding 
to 7; T q ~t >m j is the Taylor operator expanding the integrand in small masses {m^} 
and external momenta {g 7 } of the subgraph 7 (before integration); "o" inserts the 
subgraph expansion in the numerator of the integrand Fq / 7 . The sum goes over all 
subgraphs 7 which (a) contain all lines with large masses, and (b) are one-particle 
irreducible relative to lines with light masses. 
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The following types of integrals occur in the asymptotic expansion of the muon 
AMM in the Standard Model: 1. two-loop tadpole diagrams with various heavy 
masses on internal lines; 2. two-loop self-energy diagrams, involving contributions 
from fermions lighter or heavier than the muon, with the external momentum on 
the muon mass-shell; 3. two- loop self-energy diagrams with two or three muon lines 
and the external momentum on the muon mass-shell; 4. products of one-loop self- 
energy diagrams on-shell and a one-loop tadpole with a heavy mass. Almost all of 
these diagrams can be evaluated analytically using the package SHELL2 For 
our calculation we have modified this package in the following way: 1. There are 
no restrictions on the indices on the lines (powers of scalar denominators). 2. More 
recurrence relations are used, and the dependence on the space-time dimension is 
always explicitly reducible to powers of linear factors. 3. A new algorithm for sim- 
plification of this rational fractions is implemented. These modifications essentially 
reduce the execution time (in some cases, down to the order of a hundredth). 4. 
New programs for evaluating two-loop tadpole integrals with different masses are 
added. 5. New programs were written for the asymptotic expansion of one-loop 
self-energy diagrams (relevant for renormalization) in the large-mass limit. 

3. The toy-model 

As the first step, we concentrate on a "toy" model, a "slice" of the Standard Model, 
involving a light charged spinor 'J, the photon A lll and a heavy neutral scalar 
field The scalar has triple (<?) and quartic (A) self-interactions, and the Yukawa 
coupling to the spinor (y). The Lagrangian of the toy-model reads (in the Euclidean 
space-time): 

(d + mj$ + ie^A® - y$^^ (2) 

where e is the electric charge and a is a gauge fixing parameter. 
The main aims of the present investigation are the following: 

1. Verification of the consistency of the large-mass asymptotic expansion with 
the external momentum on the mass shell of a small mass. In particular, we check 
the cancelation of all threshold singularities, that appear in individual diagrams 
and manifest themselves as infrared poles in s. 

2. Estimation of the influence of a heavy neutral scalar particle on the AMM of 
the muon in the framework of the Standard Model. 

3. Verification of gauge independence: we use the covariant gauge with an arbi- 
trary parameter a. 

In the following we analyze in some detail the diagrams contributing to the AMM 
of the fermion in our toy-model and specify the renormalization procedure. Apart 
from counterterms 40 diagrams contribute to the two-loop AMM of the fermion. 
After performing the Dirac and Lorentz algebra, all diagrams can be reduced to 
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Figure 1: The prototypes and their subgraphs contributing to the large mass 
expansion. Bold, thin and dashed lines correspond to the heavy-mass, light-mass, 
and massless propagators, respectively. Dotted lines indicate the lines omitted in 
the subgraph 7. 

some set of scalar prototypes. A "prototype" defines the arrangement of massive 
and massless lines in a diagram. Individual integrals are specified by the powers 
of the denominators, called indices of the lines. From the point of view of the 
asymptotic expansion method the topology of the diagram is essential. All diagrams 
of the toy-model which contribute to the two-loop AMM can be classified in terms 
of 9 prototypes (we omit the pure QED diagrams). These prototypes and their 
corresponding subgraphs 7 involved in the asymptotic expansion, are given in Fig. 
[p. In dimensional regularization, the last subgraphs vanish in cases 1,4,7, and 8, 
owing to massless loops. 

All diagrams were generated in symbolic form by means of QGRAF i. The 
TLAMM-program 

1. Reads QGRAF output; 
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2. For each diagram, it creates a file containing the complete FORM program 
for calculating this diagram; 

3. Executes FORM; 

4. Reads FORM output, picks out the result of the calculation, and builds the 
total sum of all diagrams in a single file. 

All initial settings are defined in the "configuration file" . It contains information 
about the use of file names, identifiers of topologies, choice of momenta, and the de- 
scription of the model in terms of the notation that is some extension of QGRAF's. 
All diagrams are classified according to their prototypes. Identifiers for vertices and 
propagators and explicit Feynman rules are read from separate files and then in- 
serted into the FORM program. For each diagram, the relevant FORM subroutines 
are called. Since the number of identifiers needed for the calculation of all diagrams 
may exceed FORM capacity, the TLAMM-program keeps for each diagram only 
those needed for its calculation. Finally, the total sum of all diagrams is put into 
one file and can be processed further by FORM. There exist several options which 
allow one to process only diagrams 1 with given numbers; 2 of a given prototype; 3 
of a specified topology; and some debugging options. 

To perform the asymptotic expansion FORM-programs have been written for 
every prototype. For efficiency of the algorithm the following is essential: 

1. The result of the calculation is presented as a series in small parameters. Care 
is taken to avoid the production of unnecessary high powers in intermediate results. 

2. For the evaluation of the Feynman integrals it is necessary to reduce scalar 
products of momenta in the numerator to the square combinations present in the 
denomina tor Most efficiently, this is done by means of recurrence relations proposed 
by Tarasov E3 in the proceeding of the AIHENP-95. 

We use the on-shell renormalization scheme. The renormalization conditions 
are as follows: 1. The electric charge e is defined in terms of the nonrelativistic 
Thompson limit of the Compton scattering. 2. The physical masses are given as the 
pole of the propagators. Since we are interested here only in the AMM, we don't 
need the wave-function renormalization. All the tadpoles are canceled by mass 
counterterms. Then at the two-loop level the interaction $ 4 does not contribute to 
the AMM of the fermion. No special functions but the logarithms appear in the 
final result: 




6 Towards automatic analytic evaluating of massive Feynman diagrams 




where S 2 = ^Cl 2 (§) = 0.2604341 with Cl 2 the Clausen function, e 2 /(47r) = 
and y, g are renormalized coupling constants. We see that 

1 The two- loop contribution to the AMM is gauge independent. 

2 The additional threshold singularities arising in the asymptotic R-operation 
have canceled. 

3 With the Standard Model values it is ~ (||) 4 . 
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